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Abstract 

Relativistic elasticity on an arbitrary spacetime is formulated as a 
Lagrangian field theory which is covariant under spacetime diffeomor- 
phisms. This theory is the relativistic version of classical elasticity in 
the hyperelastic, materially frame-indifferent case and, on Minkowski 
space, reduces to the latter in the limit c —> oo. The field equations 
are cast into a first - order symmetric hyperbolic system. As a conse¬ 
quence one obtains local-in-time existence and uniqueness theorems 
under various circumstances. 
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1 Introduction 

By far the most popular relativistic matter model has been that of a perfect 
fluid. One reason for this lies in the relative simplicity of this model. Another 
reason is that astrophysical objects, which are the objects most likely to 
exhibit effects of Special or General Relativity, are usually described quite 
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well as bodies composed of perfect fluid. Elasticity, on the other hand, is 
relevant for describing the statics and dynamics of solids which we meet in 
everyday life, and for this relativistic effects can safely be ignored. But, 
whatever experimental relevance it might have (e.g. for describing neutron 
star crusts or bar detectors of gravitational waves), it is clearly of conceptual 
interest to have a consistent theory of bulk matter which is in accordance 
with the principles of relativistic physics. 

A basic consistency check on such a theory is to see whether its equations 
have a well-posed initial value problem. Christodoulou, in the book 0, has 
announced without proof a general existence theorem, covering the case of 
elasticity, using a new notion of hyperbolicity. We show that the elasticity 
equations on an arbitrary spacetime can in a very natural manner be cast 
into the form of a quasilinear, symmetric hyperbolic system 1 . 

Relativistic elasticity has been treated by many authors. The earliest 
references we are aware of are Herglotz [7], in 1911 (!) for Special Relativity 
and Nordstrom [2J in 1916, on General Relativity. A very influential paper 
has been that of Carter and Quintana j2j, of which a clear presentation has 
been given by Elders [ 3 . A very recent application of the Carter formalism 
is by Karlovini and Samuelsson cm We also would like to mention the work 
of Maugin ra- An excellent exposition can be found in the text book of 
Soper da. Recent useful references are Kijowski and Magli [T2j , which also 
contains an extensive bibliography, Christodoulou [2] and Tahvildar-Zadeh 

Ki¬ 
ln the present work we treat elasticity as a held theory derived from a 
Lagrangian. Our basic fields are maps from spacetime into a 3-dimensional 
” material manifold” B. Objects defined on B are physically interpreted as 
properties of the material prior to the action of deformations or other fields. 

I 11 Section 2 we review general properties of Lagrangian held theories 
with Lagrangians depending on the map / and its hrst derivatives and which 
are covariant under spacetime diffeomorphisms. In Section 3 we show that 
the Euler-Lagrange equations of these theories, when written as a hrst-order 
system, can in a very natural manner be rewritten as a ’’symmetric system”, 
i.e. as a quasilinear hrst - order system where the coefficient matrices of the 
derivative - terms are symmetric. 

Now recall that, in the quasilinear case, the coefficient matrices are in 
general functions of the independent and dependent variables. The former 

^or nonrelativistic elasticity such a formulation has been found by John 0. 
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are points in spacetime. The latter variables, in our case, consist of possible 
values for the field and its first derivatives. All these together form what we 
call the ’’deformation bundle” T> over spacetime. According to the general 
theory the symmetric system found in Sect.3 is symmetric hyperbolic at a 
point of T> if there is a covector at the corresponding spacetime point so 
that a certain positivity requirement is fulfilled for the coefficient matrices. 
The existence of such “hyperbolic points” of T> requires further constitutive 
hypotheses to which we turn in Sect.4. Here we define a class of theories 
for which T> has a subbundle consisting of “unstrained” or “natural” points 
which are in addition stress - free and for which the elasticity tensor has the 
standard ’’Hookean” structure. We show that, for such materials, our equa¬ 
tions are indeed symmetric hyperbolic near natural points in the deformation 
bundle. 

In Section 5 we discuss various physical sitations covered by our existence 
theorem and the interpretation of the initial data. 

Finally (in Appendix B) we treat the non-relativistic limit c —> oo of 
elasticity in Minkowski space. As to be expected, the basic structure of the 
equations remains essentially unchanged and we obtain again local existence 
and uniqueness. These equations are equivalent to those of the standard 
nonrelativistic theory in the hyperelastic, materially frame-indifferent case 
(see Ciarlet |4| and Gurtin El)- 

2 Preliminaries 

States of relativistic continua are best described by maps / from spacetime 
( M,g) to a 3-dimensional manifold B , called ’’material manifold”. In local 
coordinates we write X A = f A (x fl ). The space B is the abstract collection 
of particles making up the continuous medium. By ’’abstract” we mean that 
we view B as a space of its own right independent of ’’physical space”. There 
are several reasons for doing so, the most important one being that there is, 
in a relativistic theory, no natural notion of ” spatial location at given time”. 

The inverse images of points of B are supposed to form a timelike congru¬ 
ence of M. Thus there is, given /, a timelike vector u^, unique up to scale, 
so that 

u»dJ A = 0 , ( 2 . 1 ) 

Given a time orientation for (M, g), u M is fixed by requiring u M to be future- 
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pointing and 


( 2 . 2 ) 


g^u'V = -1 . 

The dynamics will be given by an action principle of the forme 

S[f] = [ p(f,df\g)>J-det(g) d A x (2.3) 

J M 

with associated Euler - Lagrange equations 



The Lagrangian p is thus a function on the ’’deformation bundle” T> with base 
M x B, coordinatized by (x p , X A , F B V ) where W 4 ^ is of rank three with null 
space timelike w.r. to g. The geometrical interpretation of F A M at (x p ,X A ) 
is that of a linear map from T X {M) to Tx(B). While we have no desire to 
enter formal excesses, we do from now on use notation which emphasizes 
the difference between a general fibre point ( X A , F B fl ) and a pair of values 
(f A (x),d IJ/ f B (x)) of (X a ,F b u ) along a particular map /. The four velocity 
u p for example can be viewed as a function on the deformation bundle. A 
convenient way of expressing the rank-condition and explicitly writing down 
this function is as follows: Choose a volume form VLabc on B and require 

a = F a ,F b v F c x n ABC (X) (2.5) 

to be non-zero and 

e^ Xp u vXp ( 2 . 6 ) 

to be future-pointing timelike, where is the volume form associated 

with the spacetime metric. This means that the linear maps F A ^ preserve 
orientation. Given this orientation, there is a positive function n on the 
deformation bundle so that 


u p = 


3 \n 


e <Wa p 


(2.7) 


is normalized. 

Along a map /, nu p is a vector field on M satifying the continuity equation 


V M (nu M ) = 0 


( 2 . 8 ) 
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identically. The physical meaning of n is that of a density of number of 
particles in the state given by /. 

We now come back to the Lagrangian p which is a function on the de¬ 
formation bundle with x-dependence only via the spacetime metric g, i.e. 
p = p(X,F;g). As a further ingredient we assume covariance of the La¬ 
grangian under spacetime diffeomorphisms, i.e. that, along a map /, p 
behaves as a scalar under transformations of the coordinates xA This is 
equivalent to the following condition on p: 

Ap = W) A + (2-9) 


for arbitrary vector fields £ /J '(x). D iffeomorphism invariance has several im¬ 
portant consequences. The first is that the energy momentum tensor T py 
defined by 

obeys the identity 

- Vj/TJ/ = £ A F A tl , (2.11) 

independently of the field equations. (Notice that both the left-hand side of 
EH> and £a should be viewed as functions on the ’’second jet space” with 
elements given by (x M , X A , F A U . F B \ p = F B (\ P )))- From Equ. (12. IIP and the 
regularity condition on F A it follows that the Euler-Lagrange equations are 
equivalent to the energy-momentum tensor being divergence-free. 

Secondly, writing out the Lie-derivatives in El) in terms of partial 
derivatives, we End the equivalent relation 


f a dp 
il dF A v 


2 g vX 


dp 

dg[l\ 


Equ. EH, in turn, has the corollaries 


( 2 . 12 ) 


and 


Equ. El implies that 


dp v n 

—- u = 0 

dg' w 

(2.13) 

f a n d P n 

P \FF}\ QpA x ^ ’ 

(2.14) 

Tfj, v u = pu p . 

(2.15) 
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Thus the energy momentum tensor has the form 

T„ u = pu^u u + S I1U , (2.16) 

where the stress tensor satisfies S^v? = 0. It follows that there exists 
tab = t (ab) on the deformation bundle so that 

T^v = pUf,u u + titabF a ^F b u . (2.17) 

(The factor n in front of tab is inserted for later convenience.) In other 
words, in the rest frame of matter, T /w has no mixed space-time components 
and the Lagrangian p is the energy density of the matter. The quantity tab 
is called ’’second Piola-Kirchhoff” stress tensor in the elasticity literature. 

As a special case consider a Lagrangian which just depends on n, as 
defined above. We find that 



dn n 

dcr 2 v 

(2.18) 

with 

hfiv 9/j.v T U/jUv . 

(2.19) 

Thus 

T fu/ — pu^Ui/ T ph^ 

(2.20) 

with 

dp 

V = n ^~ p - 

(2.21) 

Since, bv (12.111). the field equations on f are equivalent to 

= 0, we 


have recovered the standard formulation of perfect fluids: one simply con¬ 
siders Vi/Tj/ = 0 as the equations of motion for u^ as a vector field on M, 
supplements them by the continuity equation (EBD and completely forgets 
about / : M —► B. 

3 The symmetric system 

The equations Sa = 0 can be written as (we suppress the ^-dependence) 

df)d,d u f B (x ) = G A (f, df) , (3.1) 


6 




where 


(3.2) 


M^ ab := 


d 2 p 


dF B u dF A , 


= M^ ba 


The quantities M^ab can be viewed as a quadratic form on the space of 
matrices m A ^ at each point of the deformation bundle. Alternatively they 
can be viewed as a map sending matrices m A fl to matrices n B u . This map, 
the ’’Legendre map”, is in what follows required to be non-degenerate, i.e. 

AB m B v — 0 implies m A M = 0. Another restriction on M^ab follows 
from Equ. (irm Namely, differentiating EH w.r. to F A p and using the 
rank-condition one easily finds that 


u^M^ab = 0 . 


(3,3) 


We now have the properties of M^ w A b necessary for rewriting m as an 
equivalent first-order system. Let us define 

W^ab W ■= u»M Xv ab ~ 2u [x M^ ba 


= u»M X u AB + u v M XiL BA - u x M ^ 


BA 


By (13.21 1 


W"". b ( a) = W^ba'*' 


(3,4) 


(3.5) 


Furthermore, since the last two terms in EH are antisymmetric in A and v, 
from m there follows that 


W^AB (X \f,df)d x dJ B = u»(f,df)G A (Ldf) 

We now replace EH) by the following first-order system: 

W^ AB W(X, F)d x F B u = u»(X, F)G a (X, F) 

-u x {X,F)d x X A = 0 . 

We can write (inxHD in matrix form as 


(W^ab {X) 0 




0 


-5 a b u 


A I 


F 

X 


B 


B 


and see that we have a symmetric system because of (13.()K 


(3.6) 

(3.7) 

(3.8) 

(3.9) 


2 The quantity W^ab^’ a l so appears in 3), but is used there for a different purpose. 
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We have just seen that, taking a solutions X A = f A (x ) of m and setting 
F A fl = d^f A (x) we obtain a solution of (K-S.7lh.8l) . The converse is given by 
the following 

Theorem 1: Suppose we have a solution pair (X A (x), F a ^{x)) of Equ.’s 
(EHTTHD and a hyper surface E, on which 
(i) d % X A has rank three 3 

(h) (aw A -^)| s = o 

Then f A = X A (x) satishes d^f A = F A il and Equ.’s (13.11) . 

Proof: Contraction of E3> with u M and O imply 

M Xu AB d\F B v = G a (3.10) 

It remains to show that F A ^ = d fi f A . Inserting EH back into El) gives 

u x M^ BA d [x F B u] = 0 . (3.11) 

Using the regularity of the Legendre map this leads to 

u x 3 [x F b v] = 0 . (3.12) 

We now write X = (A^ A ) as a triple of 0-forms and F = F A ^dx ,x as a triple 
of 1-forms. Consider the expression F — dX. From the definition of and 
ESD it follows that 

i u (F-dX) = 0 (3.13) 

Thus 

jC u (F — dX) = i u d(F — dX) = i u dF = 0 , (3-14) 

where we use (13.121) in the last equality. The condition (ii) states that the 
pull-back of F — dX to E is zero. Thus, since rC is transversal to E, F — dX 
vanishes identically, and the proof is complete. 

Before commenting on the relation between the first and second-order 
systems we introduce a further piece of notation. The rank-condition on 
F A ^ means that it has an inverse on the space orthogonal to iC. Thus there 
exists F a ^ such that 


F a ^F b „. = 5 


zB 
A > 


f\f a v = hi 


where hjj is defined by E1H. 


(3.15) 


3 Condition (i) is the same as saying that u^(X,F;g) is transversal to E 










Defining 


(3.16) 


H ab = F A fl F B u g flu . 
we can write F% down explicitly as 

F a » = H ab F b * = H AB g^F B v , (3.17) 

where H AB is defined as the inverse of H AB , i.e. H AB H B( = 

With these definitions the relation EH implies the existence of quantities 
Hab = H(ab) and Ucadb = Udbca such that 

M^ ab = ~ Habu V + U ACB dF c ^F Dv (3.18) 

Thus the regularity conditions is equivalent to both hab and Ucadb being 
regular in that Hab<a b = 0 =>■ a A = 0 and U AB cd(a CD = 0 =>■ a AB = 0. 
Next, using EH and EH, it follows that 

W lw AB (x '>u x = habuW + U ACBD F Cil F Dv . (3.19) 


Note that, compared to M ,w ABj only the sign of the first term has changed. 
The regularity of AB is thus equivalent to the requirement that u M be a 
non-characteristic covector for the first-order system EZEH- 

Now recall the notion of characteristic covectors for the second-order 
system (EH). Namely k tl is characteristic iff 

Aab = M^ABkpK (3.20) 


is singular, i.e. A = det( A ab ) = 0. If is characteristic in this sense and 
a A an associated eigenvector, i.e. A AB a B = 0 , then W^ ab^ x ' k\k u a B = 0 
i.e. k^ is also characteristic for the system (13.713.81) . 

The first-order system (13.713.81) is called symmetric hyperbolic if, in ad¬ 
dition to the symmetry of the coefficients, there exists a ’’subcharacteristic” 
covector (hypersurface element), that-is-to-say a covector k\ for which the 
matrix 

0 , (321) 


AB^k\ 
0 


-H 


AB 


U X k\ 


is positive definite in the variables ( m A fJ ,,l A ). In fact, using the variables 
m A ^u^ = — a A , m A ^ L F B,J ‘ = a AB , l A , the quadratic form corresponding to 
Em for u\ = k\, is given by 


HABOt A Q[ B + U A cBDOt Ac a BD + H AB l A l B . (3.22) 

i H A B describes the distance of particles of B in spacetime. 
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Therefore, u x is subcharacteristic for (13.713. 8 j) iff the first two forms in (18.221) 
are positive definite. Notice that these conditions are also sufficient for reg¬ 
ularity. The validity of these conditions is studied in the next section. 


4 Hyperbolicity 

We first return to the covariance condition Equ. dZH2D and claim that it is 
equivalent to the requirement that p(X, F; g ) is a function just of (X, H ) 
with H ab given by (ism To prove this assertion one first observes that 
Eau. (12. 181) implies 

p(X A , F B fl ; g vX ) = p(X A , F B „ g vX + su u u x ) (4.1) 

for all real numbers s. Setting s = 1 in Equ. ED and using that = 

F a ,F b u H ab we infer that there is a function cr so that 

p(X A , F b g uX ) = a(X A , F B „ H CD ) . (4.2) 

Using ED again in Equ. E2D we see that o has to be independent of F A ^ ll 
thus proving our assertion. 

We now write p as 

P = ne , (4.3) 

where e is a positive function with 

e = e(X A , H bc ) . (4.4) 


The quantity e is the relativistic version of the “stored-energy function” of 
standard elasticity. It is possible to factorize p in this way since the number 
density n is uniquely given in terms of ( X A , H BC ) by virtue of 

6 n 2 = H AA 'H BB 'H cc 'n AB c^A'B'C' , n > 0 . (4.5) 


Just as a pair (X A , F B U ), along a map / : M —> B, measures deformation, 
the pairs (X A , H BC ) measure ’’strain”. Thus triples (x,X, H ) might be con¬ 
sidered as forming the ’’strain bundle” over M x B in which the deformation 
bundle V is embedded via Equ. EHD- Using the definition of F A M we easily 


find that 


dF A ^ 

dF\, 


-F/F b » - H ab uFu v 


(4.6) 
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and, from the definition of u^, that 

d 


-F a *u v . 


(4.7) 


dF A „ 

From (HI we deduce that 


dn 

dF\ 


nF A » 


(4.8) 


which, using (14.614.7ft . leads to 


d 2 n 

8F B v dF A „ 


2nF A ^F B u] - uHabu^vF . 


(4.9) 


The quantities /j,ab and Uacbd entering the expression (13.181) for M ,lu ab 
using dura can now be written as 


Hab = n(eH AB + t A b) (4.10) 

di~ 

Uacbd = ^{ t abHcdFtacjFbdFtbdH AcF'2^—p^+2eF[ a[cjHb}b) , (4.11) 

where the second Piola-Kirchhoff stress tensor tab entering (Era is given 
by 

tab = 2g^g. (4.12) 

Note that hab and Uacbd are functions solely on the strain bundle. Further¬ 
more the dependence on the volume form VLabc is only via the multiplicative 
factor n. It will follow that our results in this paper do not depend on the 
choice of VLabc- 

O ° ° 

Of particular interest are points (x,X, H) of the strain bundle which are 
stress-free i.e. for which t A b vanishes. Then the regularity condition on /. i A b 
is clearly satished since e > 0, but that on Uacbd becomes delicate: Since 
the fourth term in (Era, in addition to Uacbd = Ucadb, is symmetric in 
both ( BD) and (AC), this term annihilates elements a AB with a AB = a^ AB ^ 
which the last term simply multiplies by a positive constant. On the other 
hand the second term annihilates elements a AB with a AB = a\ AB \ Thus reg¬ 
ularity requires a balance between these two terms, which is not necessarily 
satished in the cases one wants to consider. But there is a possible cure: One 
can try to add further terms proportional to H C [aHb\d to Uqabd so that 
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regularity is satisfied. By that same token one might hope to also satisfy 
part 2 of the hyperbolicity requirement, namely the existence of a timclikc 
covector k lt . We would, by this manouvre, change M' 1U a r into M ,J,V ar and 
accordingly W IJ ' V ar'U into W^ab^- The second-order system, however, 
would be unchanged. We would thus be able to use W^ab^ instead of 
W^ u ab^ in Theorem 1. Consequently the statements at the end of Section 
3 can be generalized by saying: The system (13.713.31) with W^ab {X) is sym¬ 
metric hyperbolic if there exist quantities A acbd on the strain bundle with 
A acbd = A cadb, but also A acbd = ~^adbc, in such a way that 

Uacbd = Uacbd + A acbd (4-13) 

is positive definite. A necessary condition for this to be the case is that 
Uacbd is positive definite on rank-one elements, i.e. on elements a AC of the 
form 

a AC = A A n c , (4.14) 

since that property remains unchanged under mm - Thus the question 
is whether rank-one positivity, called Legendre-Hadamard condition in the 
time independent theory, is sufficient for the existence of a suitable A abcd- 
The answer, in three or more space dimensions, is in general ”no” (see Ball (T] 
and references therein). Luckily, in the situation we shall presently consider, 
the answer is affirmative 5 . 

We firstly imagine there to be given a positive definite metric G AB {X) 
(usually taken flat in standard elasticity) on B playing the role of “zero 

o ° ° 

strain”. Points (x,X,F) on V will be called “natural” or strain-free when 

F A fl F B u g fll '(x) = G AB (X ) 6 - Next suppose that the stored energy func¬ 
tion e depends on X only via G AB and is covariant under diffeomorphisms 
of B. It is not difficult to see that this implies that e only depends on 
the principal invariants of the linear map 7i with components given by 
(H) b = H AC Gcb{X). In that case it is easy to see that M^ab transforms 
tensorially both under transformations of spacetime and body coordinates, 
and in that case the full held equations can be written AP 1 'abV u f B = 0, 

5 Christodoulou claims that his new hyperbolic theory just requires rank-one positivity 
(see 0). The standard theorems for nonrelativistic elasticity in the second-order formu¬ 
lation also only require rank-one positivity, see |S]. 

6 We note that there is for general spacetimes no cross section X A = f A (x), F A /1 (x) = 
{dfj,f A )(x) consisting of natural points. The existence of such a map / : M —> B requires 
the flow of the associated vector field on M to be Born-rigid (see H3) 
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where V M is the appropriate ’’double-covariant” derivative. (We could, with¬ 
out altering our results, add some thermodynamics by allowing e to depend 
on an additional scalar function on B , the entropy density.) 

Finally we assume that the stress tab vanishes at natural points. We now 
expand e at natural points or, equivalently, at 7i = X, where 1 is the identity 
map. It follows that there are constants m, q , r such that 

777 

e = m + —[q tr{H - if + r{tr{H - l)) 2 } + 0((H - if) (4.15) 
8 


where m > 0 is the rest mass per particle. (In the prestressed case there 
would, in Equ. TTHt i. appear a linear term of the form p(2n) 1 (tr(7i — !))■ 
If the background pressure p is positive our results below continue to hold 
but the expressions for the phase velocities change.) 

The connection between q, r and the Lame “constants” is A = m n q and 

° 7 

fi = m n r. 

From and gut we deduce 


° o 

UACBD= VI U 


qGAcGsD + ‘2i"Gd(cGa)b + 2G a[cG d]b 


(4.16) 


o 

One easily finds that U is rank-one positive iff r > 0 and 2r + q > 0. We set 



2 o 2 2 

q — c 1 — 2c 2 , r = c 2 . 

(4.17) 

Inserting (|4.1(i[) and 

tab= 0 into M^ab we find that A 

= det( A ab ) = 

ABkuku: in a frame adapted to -u M , is given by 



(m °nf{9^k,K)(g^k,k v f 

(4.18) 

where 

r,f + (l-i)>S 1 ' 

c i 

(4.19) 

and 




r = f + (i- 7 )^?. 
c 2 

(4.20) 


7 Since in our formulation the volume element fi is independent of the metric G AB , 

O 

there is in general no need for n to be constant. Thus the Lame coefficients could depend 

O 

on X , but not their quotient. The phase velocities of sound are constant, as we see below. 
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Thus kn is characteristic for the second order system at natural points if 

1 2 

either Ay is null w.r. to 9 or Ay is null w.r. to 9 fJ,u . For the associated 

vectors a A with A ab& B = 0 we get that a A is proportional to F A M A; M and 
2 12 

a A is determined by a A a b Gab — 0. Thus we have a longitudinal mode 
propagating at phase velocity c\ and two transversal modes propagating at 
phase velocity o 2 . 

We now come to the question of symmetric hyperbolicity of (I3.7I3.8[) . pos- 

O 

sibly after modifying Uacbd by changing the factor of 2 in front of Ga[cGd)b 
in (TTTBt) into a generic constant a. Setting a = 4 c\ — 25 we find the following 
Theorem 2: Let 0 < 5 < ^ and 0 < 5 < 2c|. Then the system (I3.7l3.8j) 
is symmetric hyperbolic at natural points. 

Proof: We show that Ay =u ll is timelike. Keeping in mind (I3.22J) . this is 
equivalent to the positivity of S(a,a) given by 

S(a , a) = [(c? — 2c 2 2 )GacGbd + 2 c\Gd(cGa)b (4-21) 

+(4c 2 2 - 2 5)G A [cG D]B \a AC a BD > 0 

Decomposing 

a AB = u AB + k ab + -G ab (4.22) 

3 

where iv AB = uM B \ k ab = k^ ab ^ and k ab Gab = 0, we see that 

S(a, a) = (cl - ^-)k 2 + 5k A bk ab + (2c^ - 5)5u A b^ AB , (4.23) 

o 

from which our assertion follows. 

It follows from the above theorem that the system (rTTETHP has a well- 
posed initial value problem. But the allowable Cauchy data are quite re¬ 
stricted. Namely they have to be near ones for which the matter is initially 
static. The general allowable Cauchy data can be inferred from the following 
Theorem 3: Suppose the equation of state e(H AB ) is of the form (14.151) 
with 0 < | c\ < c\. Then the system (13.713. 8 1) with 0 < 5 < 2c\ is symmetric 

O ° ° 

hyperbolic at the natural points on the deformation bundle (x,X,F) i-e. 
ones where F A ^ F B u g fJ ' u (x) = G AB (X)- Covectors Ay at x are timelike iff 

g^kyjku < 0 and u ,J k IL < 0 (4.24) 


14 












where 9^ v is given by inn . We remark that the additional hypothesis, i.e. 
that |c| < cf , is not necessary for hyperbolicity but convenient for a complete 
characterization of timelike covectors. This condition is however physically 
entirely reasonable since it is equivalent to that the ’bulk modulus” given by 
k = be positive. 

Proof: Recall the definition (Id.411 of W ,1IJ ab^ in terms of M^ab- In 
the case at hand M lw ar at natural points is given by 

— -M >,v ar = —Gab U>1 uU + [(cf — 2c^)GacGbd + 2cg G d ^Ga)b (4-25) 
m n 

+(4c£ - 2 5)G A[c G d]b \ F c » F Du 

We have to study -GW(k) = -^WA B ^k\m B u m A v . Decomposing 

mn mn 

m A ll = a A u^ + F A f,a B A (4.26) 

and 

k^ = lu u^ + F A fj,k A (4.27) 

there results 

— -W(k,k) = cjaAa A + (c^ — 5 )u>k 2 + 2 (S»cf)na A k A + 2 c\<jjoiabo^ AB ^ (4.28) 
m n 

+5ioaABCt BA + Ac\a A k B oi[ab] ~ 25a A k B ctAB ■ 

Here indices are raised and lowered with Gab■ In Appendix A it is shown 
that W is positive definite iff 

uj 2 

- k A k A > 0 , (4.29) 

Cl 

from which the Theorem follows. 

We now take up the discussion of characteristic covectors of Sect.3. It is 
easy to see from (EH) that the characteristic covectors, as defined there for 
the 2nd-order system, are the same for M^ u ab as for M^ u ab ■ It remains 
to find the characteristic covectors for the lst-order system (13.713.811 . Due 
to the block diagonal form of (13.713.811 the relevant determinant is given by 
— (rA/qJ 3 !), where D is the (12 x 12) - determinant 

D = det(-^W^AB W k x ) . (4.30) 

m n 
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The determinant, in a frame adapted to u is (a — cf, b = cf): 

/ p lu , 0 , 0 , a ki , (5 — b) Ic 2 , {5 — b) ks , b k% , (a — 6) ki , 0 , b k 3 , 0 , (a — 6) ki ^ 
0 , p w , 0 , (a — h) k 2 , b ki , 0 , (h — b) ki , a fcg, (6 — b) k 3 , 0 , b k 3 , (a — 6) k 2 
0 , 0 , pco , (a — 5) k 3 , 0 , 6 ki , 0 , (a — h) k 3 , b k 2 , (5 — b) k 1 , (5 — b) k 2 , a k 3 
akx , (a — 5) k 2 , (a — <5) , a a;, 0 , 0 , 0 , a a; — h a;, 0 , 0 , 0 , a w — § w 

(5 — b) k 2 , b Ay , 0,0, bui, 0,-6 w + h w , 0,0, 0,0,0 
(5 — b) k 3 , 0 , b Ay , 0,0, buj, 0,0,0, —6 w + h w , 0,0 
bk 2 , (S — b) ki , 0,0 ,— w + h w , 0,6a;, 0,0, 0,0,0 
(a — 5) kj , ak 2 , (a — 5) k 3 , aw — 5 uj , 0, 0, 0, aw, 0, 0, 0, a w — h w 
0 , (h — b) k 3 , b k 2 , 0, 0, 0, 0, 0, buo , 0, —& w + h w , 0 
b k 3 , 0 , (h — 6) Ay , 0 , 0 , — 6 w + h w , 0,0,0, few, 0,0 
0 , bk 3 , (5 — b)k 2 , 0, 0, 0, 0, 0, — 6 w + h w , 0, bcu, 0 
y (a — S) ki , (a — 5) k 2 , a k 3 , a lu — hw,0 ,0,0 ,aw — hw, 0,0,0, aw j 

which gives (using Maple) 

D = —w 6 h 5 (2 b—S) 3 (—2h+3a) (—w 2 p+6 Ay 2 +& k 2 2 +b k 3 2 ) 2 (— w 2 p+a k 2 +a k 2 2 +a k 3 2 ) 
or 

D = -5 5 K^) 6 (2c^ - h) 3 (3c 2 - 2 6)(g^k p k u )(g^k p k u ) 2 . (4.31) 

Note that A appears as a factor in D, as it has to be. 

Finally we have to find Cauchy data. Such are given by a smooth hy¬ 
persurface £ e M and on it data X(x), F(x) satisfying conditions (i,ii) of 
Theorem 1 so that the conormal of £ is everywhere timelike in the sense 
of the symmetric hyperbolic system m ■ The easiest way to achieve this is 

o o o 

as follows: Pick an arbitrary £ C M and V a point on £. Choose (A, F) so 
that G AB (X) =F A fl F B v g^(y) and, in addition, so that the conormal n p 
of £ at V is timelike w.r. to 9^ = + (1 —4) u 11 u for -u M given in terms 

of F A P by > = p A v p B x p c p h A BC■ 

3 in 

Note that this condition is satisfied automatically when cf < 1 and £ is 
spacelike w.r. to g, iv and can not be satisfied when cf > 1 and £ is timelike 

O 

w.r. to g pi> . In the other cases it can be satisfied by a suitable choice of F A V - 

— — ° — o 

We now choose a function / : £ —> B so that f A |o —X A and dif A |o =F A i- 
Choose F A P on £ so that F A P =F A M at £ and F A { = dJ A everywhere. Then, 
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in a sufficiently small neighbourhood 0 of !J, the field equations EB with 
initial data f A = f A and d^f A = F A ^ on E have a unique local solution. 

5 Discussion 

We want to discuss certain physical situations in which our results show 
existence and uniqueness of solutions of the elasticity equations. Let us stress 
that in the present paper we are concerned only with local questions. In our 
formulation of elasticity theory we have as yet no way to solve the boundary 
initial value problem corresponding to the motion of a finite elastic body, i.e. 
where the normal component of the stress tensor at the (free!) boundary of 
the body is required to be zero. Hence we have to consider just parts of the 
body or infinitely extended bodies. 

First we consider elasticity in Special Relativity. Hence our spacetime 
metric is Minkowski space, i.e. g /iu = rj llL/ . Furthermore we assume that 
the coordinates are inertial, then t = const is the natural Cauchy sur¬ 
face. We assume the body metric to be flat, i.e. G AB = 5 AB . The map 

o 

/: (t, x 1 , x 2 , x 3 ) i—* (A" 1 = x 1 , X 2 = x 2 ,X 3 = x 3 ) is a solution for the elastic¬ 
ity equation if the Lagrangian is given by G3SD- The interpretation of this 
solution is that of a relaxed body at rest at all times. The obvious physical 
question is to deform the body slightly at t — 0 and ask for solutions. Any 

O 

map / sufficiently near to / at t=0 defines data f A , F A M at t — 0 satisfy¬ 
ing the positivity requirement via Theorem 3 in Section 4. We thus obtain 
existence locally in time for pieces of the body or for an infinitely extended 
body. 

Consider next a general spacetime to be given and in it a spacelikc hy¬ 
persurface E. If the geometry on and near E is ’’close” to that of Minkowski 
space and the data close to the ones in the previous paragraph, we obtain 
again a solution local in time. 

In an arbitrary spacetime we can, following the procedure at the end of 
Sect 4, construct solutions locally both in space and time. 

Consider as a further example a spacetime of the form 

ds 2 = — dt 2 + gik(t, x j )dx l dx k (5.1) 

Let i/; be a diffeomorphism from t — 0 onto the body and define the body 
metric G AB such that i/j becomes an isometry. Then the system becomes 
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symmetric hyperbolic for the data f A = iJj a , d^f A = 0 at t = 0, or data 
near-by. We obtain solutions globally in space, locally in time. 

Acknowledgments: We thank John Ball for very helpful suggestions 
and Jurgen Elders for his detailed comments on the manuscript. 


A Appendix A 

We start with the expression mB and decompose a AB in the form 

AB AB | ,.AB | ^ r-iAB / \ -i \ 

a = IV + K + —Lx (A.l) 

o 

where u AB = u\ AB \ k ab = k,( ab ) and k ab Gab = 0. There results 

— 5 -W(k)(m,m) = LuaACt A + 2(5 — cl)a A kAK (A.2) 

m n 

+2ujcIgjabw AB + 4c 2 2 cuABa A k B 

h\j 

-2 5(uj A b + kab + -h A B)a A k B 

+(c\ - 5)um 2 + 5k(-uj A bu ab + k A bk ab + ^—) 

o 

We now proceed as follows: we first eliminate all linear Ica - terms by sub¬ 
stituting a a with Pa given by 


1 26 1 6 

(3 A = a A + —(— ~ c\)nk A + ~{pc 2 2 - S)co A Bk b - K A Bk 

U J UJ U! 


B 


Next we eliminate terms linear in kab by setting 
kab = Akab 


3k 2 2 J 

- j)k A k B 


6nkck c 


Aco 


25 1 
— 3 lAB 


Alu 


(2cl - 5)k( A VB)ck' 


c 


for some constant A p 0. Thirdly, note the identity 


(A.3) 


(A.4) 


2 K A BCK ABC = K A BK A *k C k U - ^K A Bk^K A C k 


AB 


C 


,B ..A j„C 


(A.5) 
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where kabc is the trace-free part of fiyqn&c]- he. 

K’ABC ■= ^A[B^C] + -hA[B^‘C]D^ D (A.6) 

We use (IA.6I) to eliminate kab^ b n A ck c in favour of the other two quantities 
in (1A.5I) . We do the analogous thing for lvab , using 

3tU[ABkc]^ AB k c ^ = ujAB^ AB kck c — uj A sk B to ac k° (A.7) 


The last two operations leave us with just having to worry about the signs 
of the terms proportional to k 2 , kab^ ab and ojab<^ AB ■ We have to choose A 
such that all the signs are positive and, for the ’’only if’ -direction, have to 
make an optimal choice in terms of the allowed range for kAk A . Our choice 
is A = ^(cf — y). We finally obtain 

1 


m n 


W(k)(m,m ) = 


2 ^ 

ujf3 A f3 A + Ml - ic^)kABk 

oCi 


AB 


(A.8) 


26' 2 lv , 1 kck c 


Cl 


3 

25ujc\ 


UJ 


-)k A bK 


AB 


2h w l k c k c %/ 3 k n k D 






■ )( 2 i + 


~ 2 


)K 


+u(2cl - 5 )D{^-^-)uabu ab 


UT 


3(2 cl 


<5) 


2u 


1 + 


he? 


2 uA(c? - f) 


CO[ABk C ]UJ [AB k C] 


4h 2 


+ ——KABCK 
3uj 


ABC 


where 


D(v) = 1 

Observe that -D(O) > 0 and 


2c| 


(2c: 


2u 2 


4 u 2 (cl 


S)5c 2 2 


(A.9) 


d A 

c 2 


5 2 


4 r 2 

3 C 2 


4c| cj 


2 S 
3 


> 0 


(A.10) 


Consequently D(rj) is positive for all G [0,4], since 4 < 4 and D(rj) is 

C i C- 1 C 2 

> 0, W(k)(m,m) 


monotonically decreasing for rj > 0. Thus for 4r 


k^k 
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is positive, except possibly when loab, kab, K and (3a are all zero. But then 
kab = 0 = ola- 

Suppose, conversely, that 4 — kA \ < 0. Then choose ujab = 0, k non¬ 
zero, kab such that Rab = 0 and a .a such that /3a = 0. Clearly kab is 
non-zero of the form ’’trace-free part of &a 4 b’\ whence the last line in (I A. 811 
is zero,as are the first and fourth line. But the second and third line are 

non-positive. Thus, COD is positive def ini te iff 4 — kA^ A > 0 . 

C 1 

B Appendix B 


Here we outline how nonrelativistic elasticity is recovered by taking the limit 
c —> oo in our equations when ( M,g) is the Minkowski space 8 . We have 


g^dx^dx 1 ' = —c 2 dt 2 + 5ikdx l dx k 

(B.l) 

g*'d li d v = -±(d t ) 2 + 6 ik d i d k 
c 

(B.2) 

and 


u»d, = - - T (d t + v%) 

(B.3) 

1 Vi 

Undx^ (—cdt H— -dx l ) 

vWf 

(B.4) 

where V{ = 5ikV k and v 2 = VjV 1 . The vector v l is determined from F A u dx M = 

F A dt + f A idx l by 


F a + F A i v i = 0 

(B.5) 

The tensor h^ u = dp + u^u u goes in the linit c —> oo to 5£ 

— where 

= d t + v l di 

(B.6) 

and 


r^dx 11 = dt 

(B.7) 

The quantity Fa^O^ tends to Fa 1 d r and 


F A »F A V = 6? - v»r u , F a ^ = 5 A 

(B.8) 

8 A related discussion can be found in T5j. 



20 








The quantity H 


AB = F A fl F B u g ,iU tends to 

K ab = F A i F B j 5 ij (B.9) 

Finally, the equation of state e(X A , H BC ) can not be expected to have a 
finite limit as c — > oo, but e(X A , H BC ) has a limit, where 

e = me 2 + e(H AB ) = me 2 + e{I< AB ) + 0(\) (B.10) 

c 2 

We omit the term me 2 ( which diverges as c —> cx)) from the Lagrangian 
density (IB.101) because it contributes nmc 2 to the Lagrangian density which 
has vanishing Lagrangian derivative and thus leaves the field equations un¬ 
changed. We insert (14.1014.111) into (EHD and take the limit c —> oo to finally 
obtain 


M^ u ab = —nmJ<AB v^v v +n 

where F Afl = F A ^h^ v with 


tabKcd + tacKbd + tbdKac + 2 


dr A c 

dK BD 


pCfi pDu 


(B.ll) 


Furthermore tab — 2 
6 n 2 - 


h^d,d v = s lk dA . 

qj^ab and 7i is defined by 
= K aa, K rr 'K cc '{I arc Ma'B'C' 


, n > 0 


(B.12) 


(B.13) 


Thus the structure of the nonrelativistic equations is very similar to that of 
the relativistic ones. In particular, we can obtain existence (locally in time) 
by essentially the method described in the body of the paper. 

In nonrelativistic elasticity it is common to use as the basic field variable 
4> l (t , X) defined by 

f A (t,<P(t, X)) = X A (B.14) 

We state without proof that the equations derived above, when rewritten 
in terms of (j), are equivalent to the ones usually considered in the frame- 
indifferent, hyperelastic case (see Gurtin [BJ). 
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